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I .  INTRODUCTION 


Since  the  early  1960's,  when  gas  lasers  with  visible  illumination 
were  first  available,  the  most  casual  observer  of  radiation  from  these 
devices  has  been  able  to  note  with  interest  the  granular  pattern  come 
to  be  known  as  a  speckle  pattern,  which  accompanies  coherent  radiation. 
In  the  past  fifteen  years,  a  large  amount  of  surprisingly  diverse  work 
has  been  undertaken  which  nevertheless  falls  under  the  heading  of 
speckle  phenomena.  The  speckle  effect  has  been  used  to  describe  spatial 
noise  in  acoustic,  microwave,  and  optical  holographic  reconstruction, 
the  diffraction  pattern  resulting  from  coherent  optical  scatter  from  a 
diffuse  surface,  the  granular  noise  present  in  atmospherically  degraded 
telescope  images  and  even  the  corresponding  irradiance  distribution 
within  the  receiver  pupil.  The  thread  common  to  each  of  these  studies 
is  that  when  radiation  of  sufficient  coherence  is  scattered  from  or 
transmitted  through  a  surface  or  volume  characterized  by  random  phase, 
a  diffraction  pattern  of  a  stochastic  nature  is  generated.  It  is  not 
necessary  for  the  source  of  radiation  itself  to  be  highly  coherent;  by 
means  of  propagation  la  the  van  Cittert-Zernike  theorem)  the  coher¬ 
ence  of  the  radiation  may  attain  levels  sufficient  for  the  formation 
of  speckle  patterns  in,  for  example,  stellar  imagery.  Indeed,  the 
history  of  speckle  patterns  is  now  a  century  old. 

In  general,  the  nature  of  a  received  speckle  pattern  may  be  due  to 
the  characteristics  of  the  radiation  source,  a  surface  from  which  the 
energy  may  have  been  scattered,  or  the  characteristic  of  the  medium 
through  which  the  radiation  has  propagated.  In  this  paper  we  wish  to 
utilize  in  a  device  some  characteristics  of  the  second  type,  the  nature 
of  the  received  speckle  statistics  as  they  relate  to  the  surface  through 
which  coherent  light  has  been  scattered.  There  are  still  two  areas  of 
concern  within  this  subset:  the  phase  characterization  of  the  scatter¬ 
ing  surface  and  the  amplitude  distribution  over  the  surface.  It  is  this 
latter  aspect  of  the  problem  that  we  are  addressing  here. 

In  a  previous  paper^,  the  subject  of  image  information  by  means  of 
speckle  pattern  processing  was  discussed  with  a  view  to  the  unifying 
aspects  of  a  series  of  different  detection  and  processing  schemes.  As 
noted,  one  speckle  processing  method  having  roots  in  the  Hanbury  Brown 
Twiss  experiment,  is  called  irradiance  correlation.  To  review  this 
technique,  we  refer  to  Figure  1.  Coherent  light  illuminates  a  surface 
having  a  random  phase  distribution.  In  the  far  field  the  resulting 
irradiance  distribution  is  detected.  We  take  two  points  in  the  detec¬ 
tion  plane,  form  the  irradiance  product  and  seek  the  expectation  value 
(ensemble  average)  of  this  product.  The  domains  available  to  us  by 
which  we  can  infer  this  average  are,  of  course,  the  time  and  space 


1.  P.H.  Deitz,  "Image  Information  by  Means  of  Speckle-Pattern  Proc¬ 
essing",  Journal  of  the  Optical  Society  pp.  279-285  (1975). 
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domains.  That  is  to  say,  if  the  received  irradiance  (speckle)  pattern 
does  not  change  in  time,  we  may  translate  about  the  pattern  forming 
products  of  irradiance  pairwise,  all  for  the  same  relative  separation 
and  orientation,  and  compute  the  average.  This  is  a  spatial  average. 

Or,  if  the  source  phase  characteristics  vary  in  time,  we  may  remain 
at  the  same  two  points  in  the  detection  piano  while  the  pattern  changes 
in  time  and  forms  an  average.  This  is  an  average  in  the  time  domain. 

Or  it  is  possible  to  form  a  combination  of  averages. 

As  discussed  previously,^  the  second  moment  in  irradiance  for  the 
experiment  described  can  be  expressed  as 

<  I(xp  I  (£2)  >()“!<  v(xp  V*(X2)  >  (  )  I" 


00 

y  1(0  exp  ^  (2il  ■  2^2^  *  il  ‘^i 
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where  v  represents  the  spatial  fourier  transform  of  the  source  electric 
field,  X  and  £  represent  positions  in  the  receiver  and  source  planes, 
respectively,  k  is  the  wavenumber  of  the  light,  and  is  the  distance 

between  the  source  and  receiver  planes.  The  parentheses,  (  ),  imply  a 
domain,  not  yet  specified,  in  which  the  averaging  takes  place.  Equation 
(1)  relates  the  second  moment  in  irradiance  to  the  spatial  power  spectrum 
of  the  source  radiance  distribution.  Since  the  power  spectrum  of  a  func¬ 
tion  does  not  uniquely  define  the  function,  even  a  noise-free  measure 
cannot  be  used  to  infer  the  radiance  distribution  of  an  arbitrary  ob¬ 
ject.  The  problem  is  further  complicated  because  of  the  relatively  poor 
signal-to-noise  characteristic  of  an  incoherent  detection  technique. 

Thus  irradiance  correlation  as  a  general  means  of  source  radiance  esti¬ 
mation  is  not  particularly  useful. 

However,  in  the  context  of  a  detection  scheme,  this  means  of  infor¬ 
mation  processing  may  have  more  utility.  Rather  than  inferring  the 
radiance  distribution  for  an  arbitrary  source,  a  set  of  orthogonal  source 
functions  is  chosen  where  the  orthogonality  is  defined  in  the  spatial 
power  spectrum  of  object  radiance.  Then  the  signal  processing  problem 
at  the  receiver  reduces  to  one  of  detection  rather  than  estimation.  A 
communication  process  is  therefore  defined,  that  might  be  termed  a  spa¬ 
tial  multiplex,  spatial  diversity  scheme.  The  spatial  multiplex  aspect 
arises  because  of  the  necessary  presence  of  a  random  phase  screen 
(ground  glass)  at  the  source  to  guarantee  the  spatial  incoherence  of  the 
radiation  eminating  from  the  object.  The  presence  of  this  element  de¬ 
stroys  the  normal  diffraction  pattern  of  the  object  in  much  the  same 
way  that  a  spread  spectrum  modulation  technique  disguises  the  temporal 
character  of  transmitted  information.  The  spatial  diversity  arises,  of 
course,  from  the  multiplicity  of  detection  points  used  in  the  receiving 
process. 
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In  most  communication  systems,  the  message  coding  is  handled  by 
temporally  modulating  the  total  intensity  of  the  transmitted  beam. 

There  are  two  principle  disadvantages  to  this  process:  (1)  the  turbu¬ 
lent  atmosphere  in  the  propagation  path  also  modulates  the  beam  in  space 
and  time,  causing  errors  in  the  detection  process  and  (2)  the  presence 
of  aerosols  along  the  optical  path  causes  some  of  the  optical  energy  to 
be  scattered  out  of  the  beam.  This  off-axis  scatter  can  be  observed  at 
large  angles  to  the  initial  direction  of  beam  propagation  and  by  monitor¬ 
ing  the  irradiance  fluctuations,  the  message  being  transmitted  can  be 
inferred. 

In  the  technique  being  discussed  here,  the  information  being  trans¬ 
mitted  is  coded  in  the  cross  sectional  (spatial)  character  of  the  optical 
beam,  not  the  total  energy  being  transmitted.  To  infer  the  message,  the 
beam  must  be  detected  at  many  points  within  its  extent  and  the  detected 
signals  then  properly  processed.  Radiation  scattered  from  the  beam  does 
not  retain  the  spatial  character  necessary  for  decoding  since  it  is  only 
related  to  the  total  energy  being  transmitted.  In  addition,  because  the 
detection  does  not  take  place  in  the  fourier  transform  plane  of  the 
receiver  optics,  the  method  is  relatively  insensitive  to  the  degrading 
effects  of  atmospheric  turbulence  by  comparison  with  standard  image¬ 
forming  methods. 


II.  THE  DEVICE 

For  a  description  of  the  operation  of  this  communications  device*", 
we  refer  to  Figure  1.  This  schematic  shows  a  simplified  diagram  of  the 
equipment  layout  used  to  test  the  basic  physical  relationships.  The 
light  source  is  a  laser.  A  coherent  source  of  light  is  critical  to 
this  method  since  this  approach  is  based  on  the  generation  of  laser 
speckle  patterns  in  the  optical  farfield.  Light  from  the  laser  falls 
directly  on  the  plane  labeled  "spatially  encoded  message."  This  message 
bit  is  composed  of  a  section  of  ground  glass  in  close  proximity  to  an 
aperture  function,  typically  an  isometric  bar  grating  of  some  spatial 
frequency  and  orientation.  Provision  is  made  for  inserting  another 
isometric  bar  grating  of  a  second  spatial  frequency  (same  orientation) . 

The  ground  glass  is  arranged  so  that  it  can  be  spun  by  means  of  a  motor. 

In  the  diagram,  a  lens  is  used  to  achieve  the  far  field  optics  condition 
normally  acquired  through  propagation  over  a  long  free-space  path.  The 
use  of  the  lens  affords  a  geometrical  compactness  for  this  bench  setup. 
Following  the  lens,  a  rotating  mirror  intercepts  the  beam,  deflecting 
it  to  a  bank  of  three  photomultipliers.  The  motion  of  this  mirror  causes 
the  laser  beam  to  be  swept  across  the  three  photomultipliers,  giving 
effectively  a  one-dimensional  spatial  scan  through  the  laser  beam.  This 
scan  is  made  perpendicular  to  the  orientation  of  the  bars  of  the  isometric 

Y.  US  Patent  No.  4,085,319;  awarded  April  18,  1978. 
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Figure  1.  Laboratorv  setun  to  demonstrate  comunication  scnei"e.  'rtical 
source  is  a  laser  wnose  output  is  directed  to  a  spatially  encoded  siessa'je 
(qenerally  an  amplitude  qratinq!  in  close  proximity  to  a  random,  onase 
screen  (around  glass).  Transform  lens  gives  the  fraunnoffer  .nactern  of 
the  spatial  message  at  a  foreshortened  range.  Multisided  rotating  mirror 
reflects  beam  to  oank  of  three  photomul tiol iers.  Insert  gives  details  of 
spatially  encoded  message.  Random  phase  screen  ?S  (ground  glass)  is 
present  in  aperture  H.  Amplitude  grating  .lA  (or  GC)  forms  soatial  mes¬ 
sage  bit  for  transmission  to  detection  plane. 


Figure  2.  Dual  channel  cross-correl.it1on  device.  Output  of  ohoto- 
multiplier  tube,  R,  is  cross-correlated  simulataneously  with  the 
outputs  of  tubes  A  .and  6.  Functions  are  estimated  by  multiolyinq 
and  integrating  in  each  channel  for  a  time  t.  At  the  end  of  the 
cycle,  a  comoarator  indicates  the  channel  with  the  higher  integrated 
power . 


bar  Kratinn.  The  I'M  tube  apertv>res  are  about  l/S  the  diameter  of  the 
avcrajie  speckle  size  present  in  the  far  field  to  prevent  aperture  aver 
as  ins. 

For  a  particular  isometric  bar  specific  spatial  frequency 

information  is  present  in  the  far-field  speckle  pattern.  I'his  informa¬ 
tion  in  the  far-field  speckle  pattern  is  mathematically  related  to  the 
fourier  transform  of  the  intensity  at  the  isometric  bar  statins.  When 
the  speckle  pattern  is  detected  by  the  PM  tubes  (a  squarins  process!  and 

t 

electronically  cross-correlated  (C  =  ^  Ij^I^ilt,  where  =  output  of 

o 

PM  tube  R,  =  output  of  PM  tube  A  and  t  =  some  predetermined  time 

interval),  the  resultins  correlation  function,  C,  is  related  to  the 

power  spectrum  of  the  intensity  on  the  orisinal  isometric  bar  statins 

(see  Fis.  2).  This  correlation  function  can  be  ma.ximized  for  a  specific 

bar  statins  by  separatins  a  pair  of  PM  tubes  by  a  distance  correspond  ins 

to  the  far  field  spatial  las  fbe  fundamental  spatial  frequency  of  the 

bar  statins.  If  PM  tubes  R  and  A  are  separated  properly  for  bar  statins 

"a",  and  PM  tubes  R  and  B  are  separated  properly  for  a  second  bar  statiits 

"b",  then  a  binary  communication  system  has  been  established;  while  the 

speckle  pattern  is  interrosated  for  spatial  frequency  components,  only 

one  or  the  other  is  present  in  this  particular  setup.  Thus,  if  bar 

sratins  "a"  is  present  in  the  beam,  then  the  cross  correlation  function 

produced  from  PM's  R  and  A  (C  )  will  be  a  roiximvim,  while  the  cross  cor- 

A 

relation  function  produced  from  PM's  R  and  B  will  be  a  minimum 

(noise  level).  Replacins  bar  statins  "a"  with  "b"  reverses  the  results; 

C  becomes  sreater  than  C  .  Speckle  information  is  made  available  to 

the  PM  tubes  for  detection  and  processins  by  three  means:  spatial  scan, 
temporal  scan,  and  a  combination  of  spatial/temporal  scans. 

In  the  spatial  scan  mentioned  above,  the  stound  laser 

transmitter  is  fixed,  and  the  rotatins  mirror  causes  the  nonchansins 
laser  speckle  pattern  to  move  across  the  PM  tubes  thus  civinj;,  in  this 
instance,  a  one-dimensional  spatial  scan  across  the  sjieckle  pattern. 

In  the  temporal  case,  the  jiround  ylass  (located  close  to  the  bar 
RratinR)  is  rotated,  causing  the  speckle  pattern  itself  to  change  with 
time  while  the  rotation  mirror  is  held  fixed.  Phus,  the  three  PM  tubes 
detect  at  only  three  specific  places  in  the  laser  beam,  but  .»re  provided 
with  a  changing  speckle  pattern  produced  by  the  rotating  grovind  glass. 
(Note,  however,  that  the  speckle  pattern  at  0  radians  is  reproduced 
periodically  at  2it  radian  rotations.) 

In  the  temporal/spatial  combination  the  ground  glass  and  the  mirror 
are  both  rotated,  providing  greatly  increased  statistical  samplittg.  I'he 
combined  method  has  the  advantage  of  making  the  maximum  number  of  speckle 
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statistics  available  for  detection  and  processing  and,  hence,  increasing 
the  system  signal-to-noise  ratio.  In  addition,  by  scanning  through  the 
beam,  spatial  noise  due,  for  example,  to  atmospheric  turbulence  (which 
remains  essentially  unchanged  during  the  transnus:'i''n  of  a  message  bit) 
tends  to  average  out  spatially. 

Thus,  in  this  simple  laboratory  demonstration,  one  of  two  bar  grat¬ 
ings  is  placed  in  close  proximity  to  a  section  of  ground  glass  at  the 
transmitter  and  illuminated  by  coherent  light.  If  three  photomultiplier 
tubes  are  used  in  conjunction  with  a  rotating  mirror  configuration  at  the 
receiver,  the  net  effect  is  to  detect  a  speckle  pattern  at  many  different 
points.  The  processing  electronics  are  comprised  of  a  dual -channel 
correlator  which  simultaneously  computes  a  time-averaged  correlation  of 
the  photomultiplier  outputs  of  R  with  A  and  R  with  B.  If  the  mirror  is 
rotated,  this  time  record  also  corresponds  to  a  space  average  over  the 
receiver  plane.  At  the  end  of  the  cominitation  time,  the  correlation 
values  are  compared.  Since  the  interdetector  spacings,  a  and  b,  are  so 
chosen  to  infer  the  fundamental  spatial  harmonics  of  bar  gratings  "a"  and 
"b",  a  significant  difference  in  power  at  the  output  of  the  dual  inte¬ 
grators  implies  the  presence  of  a  particular  grating. 

figure  3  shows  a  sample  of  osc i 1 loscojic  traces  giving  the  results 
of  the  dual  integrations.  In  Figure  .3a,  ruling  "a"  is  in  place  showing 
a  repetitive  accumulation  of  power  in  integrator  a.  It  can  be  noted 
that  integrator  b  is  fluctuating  about  some  background  noise  level.  In 
Figure  .3b,  ruling  "b"  has  been  inserted  at  the  transmitter.  Here,  trace 
"a"  reverts  to  background  level  and  channel  b  exhibits  large  deflections 
at  the  completion  of  the  timing  cycle.  Subseipient  circuitr>’  compares 
the  deflection  in  the  two  channels  and  signals  the  larger  of  the  two. 


Figure  3.  Oscilloscope  tracings  of  dual  channel  correlator,  a)  Grating 
a  is  in  transmitter  aperture  resulting  in  large  integrated  cower  in  chan¬ 
nel  a.  b)  Grating  b  is  in  transmitter  aperture  so  tti.t  cr«.oter  cower  is 
measured  in  channel  b. 
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From  the  above  description,  it  should  be  obvious  that  the  device 
description  shows  only  a  few  of  many  ways  in  which  the  basic  principle 
of  this  technique  can  be  exploited.  There  are  two  central  ideas  in  this 
communication  process.  First,  the  process  is  a  spatial  multiplex  one 
since  the  spatial  frequency  bit  (isometric  bar  grating,  etc.)  is  put  in 
close  proximity  to  a  rough  random  phase  screen  (ground  glass)  and  then 
illuminated  by  coherent  light.  The  effect  of  the  random  phase  modulation 
is  to  destroy  what  might  otherwise  be  a  recognizable  far-field  diffrac¬ 
tion  pattern  of  the  object  field  distribution.  Hence,  the  information 
at  the  transmitter  is  being  spatially  modulated  by  the  random  phase 
screen.  Second,  the  process  involves  spatial  diversity  detection.  For 
the  method  to  work,  the  far-field  irradiance  distribution  must  be  detec¬ 
ted  at  points  pairwise  in  order  to  compute  a  correlation  function.  But 
in  addition,  the  process  is  statistical  and  involves  making  averages. 

If  the  transmitter  is  arranged  to  send  a  time  invariant  signal  (the 
ground  glass  is  held  static),  the  only  means  for  obtaining  an  average  is 
by  a  traversal  across  the  receiver  plane.  This  could  be  effected  by 
using  a  multipoint  sensor  such  as  a  chargecoupled  device.  As  mentioned 
earlier,  statistical  fluctuations  can  be  smoothed  by  averaging  in  time 
from  readings  at  fewer  points  if  the  source  can  be  made  to  exhibit  a 
time  history  of  its  own.  This  can  be  accomplished  by  rotation  of  the 
ground  glass.  If  possible,  it  is  desirable  (for  reasons  which  follow) 
to  eliminate  the  ground  glass  rotation  from  communication  systems  in 
certain  applications. 

There  are  a  number  of  important  implications  of  this  process  which 
we  will  list: 

1.  Although  the  above  experiment  was  accomplished  by  modulating 
the  whole  transmitting  aperture  with  just  one  isometric  bar  grating, 
hence  transmitting  only  one  bit,  it  is  obvious  that  multiple  bits  can  be 
sent  by  a  number  of  schemes  either  by  segmenting  the  aperture  and  assign¬ 
ing  individual  portions  to  specific  spatial  frequencies  at  particular 
orientations  or  by  overlaying  various  spatial  signals  at  different 
angular  orientations. 

2.  Because  the  detection  is  made  in  the  receiver  plane  without  any 
image-forming  lenses,  random  phase  modulation  induced  by  an  interlying 
medium  (turbulent  atmosphere)  has  far  less  effect  in  degrading  the 
inferred  image  information. 

3.  Because  the  process  depends  upon  detection  of  speckle  statis¬ 
tics  at  many  different  points  within  the  beam  in  order  to  make  a  useful 
estimate,  the  total  energy  being  transmitted  has  no  hearing  on  the 
message  bit  (or  bits)  being  sent.  Hence,  aerosol  scattering  along  the 
path  (which  removes  energy  at  wide  angles  to  the  axis  of  the  system), 
does  not  provide  an  off-axis  observer  with  information  useful  in  in^'er- 
ring  the  information  being  sent  even  if  the  observer  were  to  know  the 
method  of  coding  completely. 


13 


III.  SLVPLEMENTAL  EXPERIMENTS 


A.  Various  Averages, 

In  order  to  clarify  these  points  and  to  demonstrate  some  of  the 
basic  operating  conditions  of  this  kind  of  signal  processing,  we  describe 
a  series  of  tests  undertaken  prior  to  the  construction  of  the  communi¬ 
cation  device  itself.  In  an  early  test,  the  configuration  illustrated 
in  Figure  1  was  utilized  with  the  following  changes.  As  usual,  an 
isometric  bar  grating  was  placed  at  the  transmitter  in  close  proximity 
to  the  ground  glass  and  held  static.  However,  in  the  detection  plane 
only  one  photomultiplier  was  used.  The  output  of  this  device  was  sent 
to  a  signal  processor  to  be  autocorrelated.  The  multisided  mirror  served 
to  scan  a  line  through  the  time-invariant  far-field  speckle  pattern. 

The  output  of  the  autocorrelator,  shown  in  Figure  4,  gives  the  autocor¬ 
relation  function  for  this  one  scan  as  a  function  of  t,  the  effective 
spatial  lag  in  the  detector  plane.  Since  this  function  is  an  estimate 
of  the  spatial  radiance  power  spectrum  of  the  object  (bar  grating),  the 
value  at  t  equal  to  zero  relates  to  the  average  power  of  the  bar  grating 
(of  no  particular  importance  to  us  here).  However,  the  harmonic  shown 


1  V  T  tu  a  ~ 

8£9  17.78  26S7 


Figure  4.  Correlation  of  one-dimensional  speckle  scan  vs.  spatial 
lag.  Fixed  ground  glass  (time  invariant  pattern).  Averaging  strictly 
spatial. 


at  a  lag  of  8.89  mm  gives  a  measure  of  the  power  present  in  the  detection 
plane  due  to  the  52  1/cm  fundamental  harmonic  at  the  tran.smitter .  If 
the  bar  grating  were  replaced  with  one  of  another  spacing,  this  harmonic 
would  have  appeared  at  another  lag,  indicating  power  at  some  other 
frequency.  Thus,  autocorrelation  measurements  in  the  detection  plane, 
which  depend  on  spatial  lags  of  a  particular  magnitude  and  orientation, 
result  in  an  estimate  of  the  spatial  harmonic  composition  of  a  given 
object  radiance  distribution. 
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The  ((oodness  of  the  estimate  is  related  directly  to  the  number  of 
independent  statistical  measures  extracted  at  the  detection  plane.  If 
instead  of  making  one  linear  scan  through  the  speckle  pattern,  two  de¬ 
tectors  are  set  at  a  given  spacing  and  the  statistical  variations  are 
introduced  by  spinning  the  ground  glass  at  the  transmitter  as  described 
earlier,  a  more  complete  statistical  description  is  found  for  this  exper¬ 
imental  configurat ion,  as  can  be  seen  by  the  results  shown  in  Figure  5. 
Here,  two  photomultipliers  are  used  in  the  detector  plane  while  the 
ground  glass  is  spun,  and  the  cross-correlation  function  is  examined  as 
a  function  of  detector  separation.  A  more  reliable  estimate  of  the 
object  spatial  spectrum  is  derived  at  values  of  the  higher  harmonics 
(marked  at  spatial  lags  of  17.78  and  26.67  mm).  However,  if  both  methods 
of  statistical  evolution  (time  change  through  ground  glass  movement  and 
linear  spatial  scanning  by  mirror  rotation)  are  utilized,  an  even  larger 
statistical  sample  is  brought  to  the  electronic  processor,  resulting  in 
the  combined  temporal-spatial  average  shown  in  Figure  6. 


Flijure  5.  Cross  correlition  of  speckle  pattern  from  two  points  vs. 
spatial  lag  (separation).  Ground  glass  in  transmitter  aperture  is 
moving,  giving  a  time-variant  speckle  pattern.  Averaging  is  strictly 
tempora 1 . 


Figure  6.  Cross  correlation  of  speckle  pattern  using  moving  ground 
glass  in  the  transmitter  aperture  and  rotating  (multisided)  mirror 
to  scan  spatially  as  well.  Averaging  is  a  combined  temporal-spatial 
process. 
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The  operation  of  the  communication  device  discussed  earlier  can  be 
more  easily  understood  now.  Since  two  bar  Rratinns  of  different  fre¬ 
quency  can  be  uniquely  identified  by  their  fundamental  spatial  periods, 
comparison  of  power  measurements  in  the  detector  plane  at  the  two  appro¬ 
priate  lags  can  be  used  to  infer  the  presence  of  a  particular  grating. 

B.  Simultaneous  Bit  Encoding. 

Since  the  direction  as  well  as  the  magnitude  of  the  correlation  lag 
is  important,  it  is  obvious  that  a  number  of  objects  can  be  used  simul¬ 
taneously  rather  than  just  one  at  a  time.  To  demonstrate  simultaneous 
bit  transmission,  three  different  bar  gratings  were  overlaid  on  each 
other,  two  at  90-degree  orientation,  the  third  at  45  degrees.  Tlie 
same  speckle  pattern  at  the  detection  plane  was  then  scanned  and  processed 
at  the  three  appropriate  orientations  by  the  spatial  average  method 
demonstrated  earlier  in  Figure  4.  The  results  of  this  trial  are  shown 
in  Figure  7.  In  each  autocorrelation  function,  the  position  of  the  first 
harmonic, together  with  knowledge  of  the  scan  di rect ion. uniquely  defines 
the  corresponding  bar  grating  and  orientation  in  the  object  plane. 
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Figure  7.  Three  correlation  functions  (vs. 
spatial  lag)  taken  for  three  cuts  through 
the  speckle  pattern,  Three  .implitude  grat¬ 
ings  are  present  In  the  transmitter  aper¬ 
ture:  a)  grating  A  (78  Hnes/cm)  is  hori¬ 
zontal,  b)  grating  B  (39  lines/cm)  ts  at 
45°,  and  c)  grating  C  (118  11nes/cm;  is 
vertical.  Distance  of  first  harmonic  from 
the  origin  Is  proportional  to  spatial  fre- 
guency  of  grating. 


C.  Sensitivity  to  Angular  Orientation. 

The  sensitivity  of  this  process  to  angular  misalignment  will  now  be 
demonstrated.  The  first  harmonic  peaks  for  each  of  the  three  scans  shown 
in  Figure  7  represent  concentrations  of  energy  a  distance  from  the  origin 
of  the  power  spectrum  proportional  to  the  spatial  frequency  of  a  par¬ 
ticular  grating.  That  is  to  say,  harmonic  B  is  found  closest  to  the 
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origin  (has  lowest  spatial  frequency),  harmonic  A  is  next  in  distance 
from  the  origin,  and  harmonic  C  is  farthest  from  the  origin  (has  the 
highest  spatial  frequency).  Since  the  spatial  extent  (i.e.,  uncertainty 
in  frequency  space)  of  the  three  first  harmonics  is  essentially  the 
same,  the  sensitivity  of  the  apparent  power  in  these  harmonics  should  be 
proportional  to  the  distance  of  the  harmonic  from  the  origin.  This 
should  be  true  since  the  absolute  traversal  in  frequency  space  for  a 
given  angular  rotation  gets  larger  as  the  frequency  (distance  from  the 
origin)  gets  greater. 

The  effect  is  shown  in  Figure  8,  in  which  the  power  measured  at  the 
first  harmonic  spatial  frequency  is  plotted  against  angular  scan  mis¬ 
alignment  (zero  degrees  corresponds  to  optimum  alignment)  in  degrees. 

The  three  curves  in  decreasing  width  give  the  results  for  the  three  spa¬ 
tial  gratings  in  ascending  frequency.  Hence  the  sensitivity  of  the  spa¬ 
tial  scan  alignment  in  the  detection  plane  with  respect  to  the  frequency 
grating  in  the  transmission  plane  is  proportional  to  the  spatial  fre¬ 
quency  being  utilized. 


OEORCCt  —  » 

Figure  8.  Correlation  magnitude  of  first  harmonic  vs.  angular  mis¬ 
alignment  in  detection  plane  ^or  three  amplitude  gratings.  In  eacn 
case  the  grating  is  oriented  vertically  in  the  transmitter  aperture. 

Zero  degrees  orientation  indicates  a  scan  in  the  detection  plane  hori¬ 
zontally.  Since  the  spatial  extent  of  the  first  harmonic  is  essentially 
the  same  for  all  the  gratings,  the  sensitivity  to  angular  alignment  of 
the  detection  plane  is  proportional  to  the  distance  of  the  harmonic  from 
the  origin.  Hence  detections  with  the  highest  frequency  gratings  are 
most  sensitive  to  angular  alignment.  — shows  data  from  39  lints.' 

cm  gntlng,  ■  from  78  lines/cm  grating,  and  —  —  —  —  from  118 

lines/cm  grating. 


D.  Insensitivity  to  Random  Phase  Modulation 


To  illustrate  another  of  our  assertions,  that  this  method  of  spa¬ 
tial  information  transmission  is  relatively  insensitive  to  random  phase 
modulation  along  the  transmission  path,  we  performed  the  following  test 
for  comparison  purposes.  Referring  to  Figure  1,  a  bar  grating  was 
placed  in  the  object  plane  and  the  ground  glass  held  static.  Immediately 
before  the  detection  plane,  an  image- forming  lens  was  placed  so  that  an 
extremely  small  portion  of  the  grating  image  illuminated  a  pin-hole  aper¬ 
ture  of  a  photomultiplier.  When  the  four-sided  mirror  was  rotated,  the 
image  could  be  scanned.  This  function  is  shown  in  Figure  9a,  above  the 
graphed  function.  Next  the  power  spectrum  of  this  image  scan  was  com¬ 
puted  and  is  shown  directly  below.  It  should  be  noted  that  this  is  the 
identical  harmonic  content  that  is  inferred  by  the  method  of  speckle 
pattern  processing.  The  presence  of  the  first  harmonic  is  clearly  seen 
in  this  function.  Next,  a  section  of  glass  exhibiting  weak  random  phase 
distortion  was  inserted  in  the  optical  transmission  link.  The  effect  of 
this  glass  was  to  distort  the  image.  This  degraded  image  was  then 
scanned  and  the  corresponding  spatial  power  spectrum  computed.  These 
results  are  shown  in  Figure  9b.  The  image  is  seriously  degraded,  so 
much  so  that  the  power  spectrum  no  longer  defines  the  harmonic  content 
of  the  object  spectrum. 
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Figure  9.  a)  Insert  shows  the  image 
of  a  diffraction  grating  using  a  lens. 
Graph  directly  below  shows  the  power 
spectrum  derived  from  standard  proc¬ 
essing  of  the  image.  Strong  harmortc 
content  can  be  observed,  bi  Insert 
shows  same  diffraction  grating  image 
when  a  weak  random  phase  screen  is 
present  in  the  optical  path.  Power 
spectrum  directly  below  shows  har- 


b) 
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To  test  the  sensitivity  of  the  speckle  method  to  this  kind  of  dis¬ 
tortion,  the  spatial  average  method  (illustrated  in  Figure  4)  was  uti¬ 
lized  in  the  following  way.  As  shown  in  Figure  4,  without  the  distort¬ 
ing  phase  screen,  the  detector  lag  is  8.89  mm  to  measure  the  power 
associated  with  the  first  harmonic.  The  strength  of  this  unperturbed 
harmonic  was  plotted  in  Figure  10  at  the  abscissa  position  0,  and  is 
labeled  as  such.  Then,  the  random  phase  screen  was  inserted  in  the 
transmission  path.  The  magnitude  of  the  correlation  function  in  general 
decreased.  The  amount  of  the  lag  (effective  detector  separation)  was 
then  adjusted  so  as  to  maximize  the  magnitude  of  the  correlation  read¬ 
ing.  An  increase  in  the  separation  resulted  in  a  new  magnitude  plotted 
for  positive  t,  a  decrease  in  negative  r.  Some  forty  readings  were 
taken  for  each  of  forty  particular  positions  of  the  phase  screen.  It 
can  be  seen  that  the  effect  of  the  random  phase  screen  on  any  one  read¬ 
ing  was  to  increase  or  decrease  the  apparent  frequency  of  the  harmonic 
content  of  the  ruling.  However,  when  the  separate  readings  were  aver¬ 
aged,  the  mean  measurement  converged  very  closely  to  the  true,  undis¬ 
torted  value.  Thus,  by  this  speckle  method,  not  only  was  the  harmonic 
information  not  lost  in  a  noisy  transmission,  its  apparent  (average) 
value  was  shifted  only  slightly  in  frequency. 


Figure  10.  Relative  magnitude  (indicated  by  height)  and  shift  for 
some  forty  spatial  frequency  readings  in  which  a  random  phase  screen 
was  moved  to  a  new  position  for  each  reading.  Position  of  spatial 
harmonic  in  the  absence  of  random  phase  modulation  is  also  shown.  Av¬ 
erage  (vertical  dashed  line)  of  all  readings  converges  to  unperturbed 
value. 
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r.  Further  Kandom  Phase  Statistics 


I'iRure  11  shows  further  data  describinj;  the  effects  of  turbulence 
on  the  operation  of  a  binary  communication  litik  as  described  initially. 
Here  an  isometric  bar  uratinj*  has  been  placed  in  tlie  transmitter  aper¬ 
ture.  In  the  detection  plane,  harmonic  information  was  estimated  at 
two  spatial  la^s  appropriate  for  bar  ^rat  inf’s  of  .SJ  anil  78  lines/cm. 

Both  power  readiiif’s  were  made  for  80  individual  positions  of  the  random 
phase  screen.  As  can  be  seen  in  l-'ifjure  11,  in  all  cases  but  one,  f>reater 
power  is  measured  in  the  channel  carrying  the  spatial  content  of  the  !’2 
lines/cm  firatiiif;  (labeled  sif>nal  in  I'ifture  11)  than  in  the  channel 
measurinf;  signal  content  at  78  lines/cm  and,  hence,  simply  labeled 
background  noise. 


Fiijure  11.  Voltage  in  each  cnannel  of  binary  communication  link  vs. 

30  realizations  of  experiment  in  wnicfi  random  ohase  screen  is  moved 
to  a  new  position  for  each  reading.  Plot  labeled  signal  snows  output 
from  channel  tuned  to  the  spatial  message  in  the  transmitter  aperture. 
Noise  signal  shows  output  of  receiver  channel  not  matched  to  spatial 
nessage.  In  all  cases  but  one  (the  sscond  reading),  power  in  channel 
tuned  to  spatial  message  in  transmitting  aperture  is  greater  than  the 
noise  channel  despite  strong  phase  modulation  along  the  propagation 
path. 


IV.  SUMMARY 

In  thi.s  paper,  wo  have  attempted  to  describe  a  particular  method 
of  signal  modulation,  transmission,  and  decoding  that  we  have  chosen  to 
call  spatial-multiplex,  spat i a  1 -d i vers i ty  optical  communication.  There 
are,  of  course,  many  other  techniques  that  may  be  used  to  transmit  opti¬ 
cal  data  of  a  sjiatial  nature,  including  direct  projection  of  transparen¬ 
cies  by  optical  systems  or  by  using  image- forming  optics  at  a  reception 
plane.  In  this  method,  no  image- forming  lens  is  vised.  Coherent  liglit 
is  used  to  illuminate  a  transparency  upon  which  a  bit  (or  bits')  is 
coded  according  to  amplitude  modulation  within  a  transmitting  aperture; 
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the  bit  is  characterized  by  some  spatial  frequency  and  orientation.  In 
close  proximity  to  the  message  bit  (or  bits)  the  optical  beam  must  be 
randomized  by  transmission  through  a  phase  screen  of  fast  phase  varia¬ 
tion  (such  as  ground  glass).  At  the  path  terminus  of  the  communication 
link,  the  resulting  optical  speckle  pattern  must  be  detected  at  many 
points,  with  this  data  being  properly  cross-multiplied  and  averaged; 
that  is  to  say  a  two-dimensional  correlation  function  of  the  received 
intensity  must  be  computed  and  examined  at  spatial  lags  appropriate  for 
the  chosen  transmission  format.  The  rate  at  which  information  may  be 
transmitted  depends  upon  many  complicated  tradeoffs  including  the  number 
of  bits  being  coded  simultaneously,  the  size,  number  and  speed  of  the 
detectors,  and  the  spatial  noise  along  the  link  itself. 

Of  fundamental  importance  to  this  technique  is  the  fact  that  a 
particular  message  bit(s)  being  transmitted  is  inferred  by  second-order 
correlations  of  intensity  transverse  to  the  optical  axis  and  is,  thus, 
independent  of  total  transmitted  power.  Light  scattered  out  of  the  pri¬ 
mary  beam  by  aerosols  does  not  contain  any  information  useful  to  off- 
axis  observers  who  are  constrained  to  view  only  volume  scatter  out  of 
the  beam.  Thus,  although  the  method  may  reveal  the  presence  of  a  com¬ 
munication  link  to  an  outside  observer,  the  method  offers  a  high  degree 
of  security  since  access  to  the  necessary  statistical  measures  can  be 
gained  only  by  physical  penetration  of  the  primary  beam  itself. 

In  addition,  because  the  reception  scheme  utilizes  only  direct 
energy  (intensity)  detection,  the  method  is  relatively  insensitve  to 
random  phase  modulation  along  the  link  from,  for  example,  the  effects 
of  turbulence.  Since  the  beam  is  detected  immediately,  without  passage 
through  any  phase  adjusting  devices  (such  as  lenses),  pure  phase  dis¬ 
tortion  is  lost  immediately  in  the  square- law  detection  process.  De¬ 
pending  on  complicated  system  parameters  which  relate  to  the  overall 
signal-to-noise  ratio,  more  than  one  bit  can  be  coded  and  transmitted 
simultaneously  by  either  segmenting  the  transmitting  aperture  and  as¬ 
signing  individual  bits  to  different  portions  of  the  aperture  or  by 
overlaying  various  amplitude  gratings  at  different  angles. 

Finally  we  note  that  although  in  the  spatial  diversity  scheme  dis¬ 
cussed  above  many  different  phase-modulated  points  in  the  transmitter 
aperture  were  illuminated,  spatial-multiplexing  can  be  accomplished  by 
as  few  as  two  points  illuminated  at  any  given  time.  Although  for  this 
situation  the  far-field  pattern  is  a  particularly  simple  class  of 
speckle  patterns  (e.g.,  a  fringe  pattern),  message  bit  assignments  can 
be  made  on  the  basis  of  point  (pin-hole)  separations  and  orientations 
giving  rise  to  fringe  patterns  in  the  receiver  plane  which  can  be  in¬ 
terrogated  by  many  detection  schemes  to  infer  corresponding  spatial 
orientation  and  frequency. 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN;  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 


Name: 

Telephone  Number: 
Organization  Address: 


